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1. INTRODUCTION 



The main model responsible for the energy release in active galactic nuclei (AGNs) is 
presently known to be the electrodynamic model dating back to the paper by Blandford and 
Znajek (1977). Within this model, the energy losses of a rotating black hole are related to 
the flux of electromagnetic energy flowing along the magnetic field lines from the black hole 
surface in the direction of the jets. The question about the magnetic field structure in the 
vicinity of a black hole (which should be generated in an accretion disk) still remains an open 
one. This question becomes particularly topical in connection with the latest observations of 
the inner jet regions (see, e.g., Doeleman et al. 2012) and with the successful launch of the 
Spectrum-R (Radioastron) Space Observatory. The latter allows spatial scales comparable to 
the size of the central black hole to be resolved (Kardashev 2009). 

By now, a wide variety of magnetic field geometries near a black hole have been considered 
in the force-free approximation within an analytical approach. These include a split monopole 
field near the horizon and far from the black hole (Blandford and Znajek 1977), a parabolic 
field near the horizon and far from the black hole (Blandford and Znajek 1977; Ghosh and 
Abramowicz 1997), and a uniform magnetic field near the horizon and a split monopole field 
at large distances (Beskin et al. 1992). In all cases, the angular velocity of the plasma nF(0) 
(which is known to depend uniquely on the magnetic field geometry) near the rotation axis 
has always been exactly half the angular velocity of the black hole ilu (see Fig. 1). However, 
the latest numerical simulations show that the condition r2F(0) = 0.5 Oh can be violated. In 
particular, McKinney et al. (2012) argue that their angular velocity profile Qf(^) near the 
black hole horizon more closely corresponds to a parabolic field, for which actually drops to 
0.3 d.t 9 = n/2. In this case, however, some of the magnetic field lines should be connected 
not with the black hole horizon but with the accretion disk near the ergosphere (Punsly 2001). 

In this paper, we study an analytical model of the black hole magnetosphere based on 
a previously unconsidered magnetic field geometry: a radial field near the horizon and a 
vertical field far from the black hole. In the second section, we give a brief overview of the 
Grad-Shafranov equation and other models of the black hole magnetosphere based on it. In 
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Fig. 1: Qp/Qh versus for a split monopole field near the horizon and far from the 

black hole (dashed line), a parabolic field near the horizon and far from the black hole (dash- 
dotted line), and an uniform magnetic field near the horizon and a split monopole field at 
large distances (solid line). Here \E'j, is the total magnetic flux through the horizon. 

the third section, we consider the model itself and compare it with the results of numerical 
simulations of the black hole magnetosphere (McKinney et al. 2012). We show that the 
derived angular velocity profile can be easily explained in terms of this model. 

2. THE GRAD-SHAFRANOV EQUATION 

The Grad-Shafranov equation describes axisymmetric stationary flows in terms of ideal 
magnetohydrodynamics (MHD). This approximation is based on the assumption about a good 
conductivity of the plasma filling the magnetosphere of a compact astrophysical object. In 
the vicinity of a rotating black hole (whose metric is also axisymmetric and stationary), this 
is provided by an efficient production of electron- positron pairs (Blandford and Znajek 1977). 
This approach is convenient in that quite a few integrals of motion, i.e., quantities conserved 
along the particle trajectory, exist in the case of stationary ideal MHD. This allows the MHD 
equations to be reduced to one second-order equation for the magnetic flux function \E'(r, ^^), 
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defining the magnetic field: 



B^^-^e.. (1) 

1t:w cuj 



Here, w = y/g^ is tiie distance to tiie rotation axis. For this choice of designations, the 
function ^(r, ^) coincides with the magnetic flux passing through a circle r, 9, < ip < 2Tr 
while the function 7(r, 9) is the total current flowing through the same circle. 
In addition, the following important properties hold. 

(1) The equation V ■ B = holds automatically. As a result, the three magnetic field compo- 
nents are defined by two scalar functions, ^(r, 9) and /(r, 9). 

(2) Since the equation B • = holds automatically, the ^(r, 9) — const hues specify the 
shape of the magnetic surfaces. 

Next, using the freezing-in condition E-|-vxB/c=0 and the assumption about axisym- 
metry, we can determine the electric field as follows (for more details, see Beskin 2010): 

(2) 

where u is the Lense-Thirring angular velocity. As a result, the Maxwell equation V x E = 
leads to the relation VJlp x V^' = 0, whence it follows that 

= ^^f(*). (3) 

The function introduced in this way means the angular velocity of the particles moving 
in a plasma filled magnetosphere, while condition (3) is Ferraro isorotation law, according to 
which the angular velocity of the particles on axisymmetric magnetic surfaces must be constant 
(Ferraro 1937). Similarly, from the Maxwell equations we can deduce that V/ x V^' = and, 
consequently, 

/ = /(*). (4) 

This means that the total electric current within a magnetic flux tube is also conserved. 

It is important to emphasize that in contrast to the nonrelativistic problem, there is a 
second family of singular surfaces associated with the accreting matter in the black hole 
magnetosphere. As a result, the additional critical condition allows an additional relation 



between the current /(\&) and angular velocity f2F(^) to be determined. In the force- free 
approximation, this relation can be written as (Thorne and MacDonald 1982) 



where Tg is the black hole radius and a is the rotation parameter. Recall that the true meaning 
of Eq. (5) is the critical condition on the inner fast magnetosonic surface that coincides 
with the black hole horizon in the force-free approximation (Beskin 2010). As a result, this 
condition allows not only the longitudinal current but also the angular velocity ilF(^) 
to be determined. 



As has already been said, several analytical models of the black hole magnetosphere were 
proposed in the literature. The first of them was constructed by Blandford and Znajek (1977). 
They considered a slowly rotating black hole for which a nonrotating black hole with a split 
monopole field was chosen as the zeroth approximation. Such a geometry can be easily real- 
ized in the presence of a thin accretion disk. In this case, the flux function ^ = ^o(l ~ cos^^) 
for 6* < 7r/2 and ^ = ^'o(l + cos 6*) for 6* > 7r/2 will be an exact solution of the Grad- 
Shafranov equation for a nonrotating black hole. The same authors considered a model mag- 
netosphere with a parabolic magnetic field in the vicinity of a slowly rotating black hole. 
The shape of the field lines for 6 < 7r/2 at large distances is described by the fiux function 
^ = 'i>Q{r/rg){l — cos 9). Since ^(r, tt) ^ const for it, this implies the presence of sources or 
sinks in the volume (and not only in the gravitating center or at infinity). Such sources can 
also be realized in an accretion disk. Finally, Beskin et al. (1992) investigated the case where 
the black hole is in the center of a well-conducting disk bounded by an inner radius b. In 
this case, the magnetic field was almost uniform near the black hole and still remained a split 
monopole one at large distance (r » 6). As can be seen from Fig. 1, the angular velocity 
f2F(^) near the rotation axis is f2H/2 in all these cases. 

On the other hand, as is shown in Fig. 2, in their recent paper devoted to numerical 
simulations of the black hole magnetosphere, McKinney et al. (2012) obtained a profile of 





3. THE PLASMA ANGULAR VELOCITY PROFILE 
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Fig. 2: The plot of Qp/Qn on the horizon versus polar angle derived during the numerical 
simulations of the black hole magnetosphere (McKinney et al. 2012). The dashed and dash- 
dotted lines correspond to monopole and parabolic fields, respectively. 

the angular velocity that not only differed from Q}i/2 near the axis but also even became 
negative here. The authors conclude that the derived profile is nevertheless closer to the 
parabolic solution, especially since the outer magnetic surfaces actually have such a shape. 
However, a significant fraction of the magnetic field lines in the parabolic solution must pass 
through the accretion disk. Consequently, not the rotating black hole but the equatorial region 
of the ergosphere will be an energy source for the corresponding magnetic surfaces (such a 
model was developed, for example, in Punsly's works; see Punsly 2001). 

Below, we will show that the results obtained by McKinney et al. (2012) are in best 
agreement with the previously unconsidered model of a black hole magnetosphere with a 
(split) monopole magnetic field near the black hole horizon and a cylindrical field far from 
the black hole. In this model, in which the bulk of the magnetic flux now passes through the 
black hole horizon, not only a coUimation of the magnetic surfaces (it will be connected with 
the fairly high density of the ambient medium) but also angular velocities smaller than 
near the axis can be obtained. 

In our model, we will use the assumption that the flow is cylindrical near the rotation axis 
far from the black hole, as is obtained in numerical simulations, so that all quantities depend 
only on the cylindrical radius w. In this case, the Grad-Shafranov equation is one- dimensional 
and can be integrated. In the force-free approximation, the solution of the Grad-Shafranov 
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equation takes the form (see, e.g., Istomin and Pariev 1994) 

d 



Substituting the expression of the toroidal field in terms of the total current 
we can rewrite Eq. (6) as 

where the following notation is used: 



(6) 

-21 /wc, 
(7) 



1 

^2(*) = c'jx^—iB^fdx. 



(8) 
(9) 



As regards the "boundary condition on the horizon" (5), it can be rewritten as 



(10) 



where 



As = — sin 6' 



rl + a^ /d^'' 



27r rl + a? cos^ 9 \de J ' 



(11) 



Substituting the latter expression for 21 {^) into Eq. (7) yields a quadratic equation for 
^^f(*): 

nl{Al - Al) + 20F^^i/^I - ^H^l - ^ = 0. (12) 
Hence, the general expression for the angular velocity can be written as 



A, 



+ 



(13) 



This form stems from the fact that the relation Ai — A^ holds for the previously considered 
solutions on the rotation axis. Therefore, we tried to avoid the quantities {Ai — A3) in 
the denominators of the corresponding expressions. Let us first consider the case where the 
magnetic field is vertical and uniform far from the black hole (^ = nvn'^Bo at r » Tg) and 
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exactly radial on the horizon (^^ = — cos^^) at r = Vg). Substituting the corresponding 
flux functions into Eqs. (8), (9), and (11), we will then obtain 

Id** 



2 



27rzu dzu TT 







A.m = . (14) 

' TT r2 + a2(l-*/*,)2 ^ ' 

Consequently, — 0,^/2 on the rotation axis. On the other hand, both analytical (Beskin 
and Nokhrina 2009) and numerical (Komissarov et al. 2006; Tchekhovskoy et al. 2009; Forth 
et al. 2011) calculations show that a central core, which, as will be shown below, can change 
significantly the situation, can exist in the jet. Consider the case where the magnetic field is 
still exactly radial on the horizon and vertical far from the black hole, but now a denser core 
of radius r^ore exists near the rotation axis. As has been shown, such a core must actually 
be formed at a sufficiently low pressure of the ambient medium, with (see, e.g., Beskin and 
Nokhrina 2009) 

""^"■"^ = ^?Mo) ■ ^^^^ 

Here, A; ~ 7in, where 7in is the characteristic Lorentz factor of the particles flowing along the 
jet axis. At distances zu < Tcore from the rotation axis, we can then write 

B, = Bo-Bo^. (16) 

c 

The following flux function corresponds to this field: 



'"core 



1 „ W' 



4 



^ ^TTZu'Bo- -ttBo^. (17) 



2 r 



2 

core 



As a result, we still have ^i(*) ~ ^3(*) ~ ^/tt near the rotation axis in the first order in *. 
Now, however, A2 will be nonzero, 

^2 



core 



(18) 



In the upshot, substituting the expressions for Ai,A2 and A^ into the general formula (13), 

we have 
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If, however, we express here rcore using Eq. (15), then 

Qp(0) ^ ?^===- (20) 

1 + ^1 + 

As we see, in the presence of a dense core, the angular velocity on the jet axis 0,^(0) can be 
smaller than Oh/2. In particular, for a mildly relativistic flow, = 1, we obtain f2F(0) = 0.41 On- 
Let us now use our model to analyze the results of the numerical simulations performed 
by McKinney et al. (2012), in which a dense core also takes place. According to these results, 
the magnetic field near the black hole horizon may be considered radial with a good accuracy. 
From the plot of Bj.{rii, 6), we can then derive the plot of the magnetic flux function ^'(rH, 0) 
and subsequently the plot of A^l'if). Since the flow near the axis at large distances may be 
considered cylindrical with a good accuracy, we may set ~ Br in this region and use the 
plot of the radial magnetic field at r = SOvg. However, the magnetic field directly on the 
axis is overestimated due to the peculiarities of the numerical method. Let us now consider a 
model magnetic field 

= 1 ^ ^2/ 2 + Bi, (21) 
where Bq, Bi and Tcore a-re the parameters of the problem. Let us choose them so that, first, 
the plot of the function Bz{9) at small 9 is close to the plot of Br{SOrg, 9) from McKinney 
et al. (2012), and, second, the total magnetic flux corresponding to this B^ should coincide 
with the total magnetic flux on the horizon. The latter condition is based on the property 
of magnetic flux conservation and the fact that much of the magnetic flux emerging from the 
black hole horizon is subsequently concentrated inside the jet, i.e., near the axis. The following 
flux function corresponds to the model magnetic field (21): 

* = Trr^e^o ln(l + ^V^re) + ^B,w\ (22) 

Such a flux function does not allow the inverse dependence vd{^) and, consequently, the 
dependences ^i(^) and A2{^) to be derived analytically. However, these dependences can be 
derived numerically and the angular velocity profile VL-p{^{'Uj) can be found using Eq. (13). 

Fig. 3 presents a plot where Op/On on the black hole horizon is along the vertical axis 
and the polar angle 9 is along the horizontal axis. The range of angles was chosen from 
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Fig. 3: The plot of Qp/Qn on the horizon versus polar angle 6 derived from the analytical 
model described here (solid line) and the plot from McKinney et al. (2012) (dashed line). 

the following considerations. At very small 6, the magnetic field obtained in the numerical 
simulations diverges, most likely due to the peculiarities of the numerical method. On the 
other hand, at large 6, the assumption that the magnetic field is vertical will break down. As 
we see, the proposed model is in excellent agreement with the numerical simulations. Note 
that is negative near the axis, which is also in agreement with the work by McKinney et 
al. (2012). 

4. CONCLUSIONS 

We investigated a new analytical model of the black hole magnetosphere based on a pre- 
viously unconsidered geometry of magnetic surfaces: a radial magnetic field near the horizon 
and a vertical field far from the black hole. We showed that in the presence of a dense core 
near the jet axis, there is excellent agreement of this model with the numerical simulations. 
And this is despite the fact that the analytical calculations were performed within the sim- 
plest force-free approximation and under the assumption that the flow was axisymmetric and 
stationary, while McKinney et al. (2012) carried out their 3D numerical simulations in the 
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full MHD version by taking into account the fact that the flows under consideration were 
nonstationary. 

We emphasize that the negative values of the angular velocity are most hkely associated 
with the difficulties of the numerical procedure near the rotation axis. Therefore, actually one 
should not expect the appearance of a region with counter-rotation near the jet axis. For us 
it was important here to only show that given the magnetic field structure near the black hole 
horizon and in the jet region, the angular velocity profile obtained in a self-consistent way 
could be reproduced using a simple analytical model. 

Good agreement between the theory and numerical simulations once again shows that 
axisymmetric stationary flows, for which quite a few analytical results have been obtained 
in the last three decades, remain a good basis for analyzing the processes occurring in real 
astrophysical sources. One of such properties is that despite the turbulent nature of the flow 
in the region above the accretion disk, the flow near the rotation axis remains fairly regular. 
Therefore, there is hope that the previously formulated simple analytical asymptotics (and, in 
particular, the assertion that the magnetic field structure near the horizon should be nearly 
radial) will also be needed in the future. 

Finally, note that a parabolic field such that a significant fraction of the magnetic field 
lines cross the equator within the ergosphere would require the existence of an energy source 
directly in the accretion disk. In the steady-state problem, such a situation is unlikely to be 
possible. In our view, the fact that the 3D simulations discussed above lead to a quasiparabolic 
structure of magnetic surfaces is related to a fairly high external pressure. As a result, the jet 
radius exceeds the black hole horizon radius only by several times. In fact, however, as can 
be clearly seen from the structure of the innermost magnetospheric regions in McKinney et 
al. (2012), only a very small fraction of the magnetic field hues pass through the equator. 

ACKNOWLEDGMENTS 

We are grateful to A. Tchekhovskoy for the provided numerical simulations, A. A. Philippov 
for a helpful discussion, and the Ministry of Education and Science of the Russian Federation 

11 



for financial support (Contract No. 8525). 



REFERENCES 

V. S. Beskin, Phys. Usp. 53, 1199 (2010). 

V. S. Beskin, Ya. N. Istomin, and V. I. Pariev, Sov. Astron. 36, 642 (1992). 
V. S. Beskin and E. E. Noklirina, Mon. Not. R. Astron. Soc. 397, 1486 (2009). 
R. D. Blandford and R. L. Znajek, Mon. Not. R. Astron. Soc. 179, 433 (1977). 
S. S. Doeleman, V. L. Fisfi, D. E. Scfienck, et al.. Science 338, 355 (2012). 
V. C. A. Ferraro, Mon. Not. R. Astron. Soc. 97, 458 (1937). 
P. Gliosli and M. A. Abramowicz, Mon. Not. R. Astron. Soc. 292, 887 (1997). 
Ya. N. Istomin and V. I. Pariev, Mon. Not. R. Astron. Soc. 267, 629 (1994). 
N. S. Kardashev, Phys. Usp. 52, 1127 (2009). 

S. Komissarov, M. Barkov, N. Vlaliakis, and A. Konigl, Mon. Not. R. Astron. Soc. 380, 51 
(2006). 

J. C. McKinney, A. Tchekhovskoy, and R. D. Blanford, Mon. Not. R. Astron. Soc. 423, 2083 
(2012). 

O. Portli, Cli. Fendt, Z.Meliani, and B. Vaidya, Astropliys. J. 737, 42 (2011). 
B. Punsly, Black Hole Gravitohydromagnetics (Springer, Berlin, 2001). 
A. Tchekhovskoy, J. McKinney, and R. Narayan, Astrophys. J. 699, 1789 (2009). 
K. S. Thorne and D. MacDonald, Mon. Not. R. Astron. Soc. 198, 339 (1982). 



12 



